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Abstract 

Recently a Hamiltonian formulation for the evolution of the universe dominated by 
multiple oscillatory scalar fields was developed by the present author and was applied to 
the investigation of the evolution of cosmological perturbations on superhorizon scales in 
the case that scalar fields have incommensurable masses. 

In the present paper, the analysis is extended to the case in which the masses of scalar 
fields satisfy resonance conditions approximately In this case, the action-angle variables 
for the system can be classified into fast changing variables and slowly changing variables. 
We show that after an appropriate canonical transformation, the part of the Hamiltonian 
that depends on the fast changing angle variables can be made negligibly small, so that 
the dynamics of the system can be effectively determined by a truncated Hamiltonian that 
describes a closed dynamics of the slowly changing variables. Utilizing this formulation, 
we show that the system is unstable if this truncated Hamiltonian system has hyperbolic 
fixed point and as a consequence, the Bardeen parameter for a perturbation of the sytem 
grows. 
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1 Introduction and Summary 



The inflationary universe model is the most successful model in explaining the origin of 
the present cosmological structures such as galaxies and clusters of galaxies. In this model, 
quantum fluctuations of an inflaton field, a scalar field driving the inflationary expansion, 
provide seed perturbations which grow and form the present cosmological structures by 
gravitational instability. During the slow rolling phase of inflation, these seed perturbations 
are streched beyond the Hubble horizon and their wavelengths stay larger than the horizon 
scale until the perturbations come back inside the Hubble horizon during the Friedmann 
stage after the inflation. The amplitudes of perturbations at this second horizon crossing, 
which have a direct relevance to the CMB anisotropy observations and provide the initial 
condition for detailed astrophysical models of galaxy formation, are determined by the 
so-called Bardeen parameter. Hence, in order to obtain information on the inflationary 
stage of the universe from observations of the present universe, we have to determine the 
behavior of the Bardeen parameter of perturbations during the superhorizon stage. 

If the cosmic matter has a regular equation of state and is dominated by a single 
component, this Bardeen parameter is conserved with a good accuracy on superhorizon 
scales pO]. However, in a realistic model, during the period between the first horizon crossing 
in the inflationary regime and the socond horizon crossing in the Friedmann regime, the 
inflaton field oscillates coherently around a local minimum of a potential and their energy 
is gradually transformed into matter and radiation which constitute the present universe. 
During this reheating phase, the equation of state becomes singular periodically |6 ; ] , and 
entropy modes can be produced[7m]]- Further, in an inflationary model with a multiple- 
component inflaton field, isocurvature modes can appear |TTJ Therefore, the Bardeen 
parameter may not be converved even approximately in realistic models, and a detailed 
analysis of its behavior is mandatory. 

On superhorizon scales, this problem can be reduced to the analysis of a spatially homo- 
geneous model for the following reasons. First, the evolution of cosmological perturbations 
on superhorizon scales is well described by that in the long wavelength limit with a good 
accuracy^]- Further, the long wavelength limit of a solution to the perturbation equation 
can be easily constructed from a homogeneous perturbation of the background universe 
model[HJ IH El ]. Hence, we only have to analyse a spatially homogeneous system, which 
is much simpler than the analysis of the exact perturbation equation with a finite wave 
number. If the background system is exactly solvable, this reduction solves the problem 
completely. However, in the case of multiple scalar fields in an expanding universe, it is 
not the case. Further, the system exhibits quite complicated oscillatory behavior. 

There is, however, one powerful method to treat such an oscillatory system. It is 
the Hamiltonian formulation in terms of the action-angle variables [Tol ]. In fact, in a 
previous paper by the present author [HI], we have developed a Hamiltonian formulation 
for a universe model dominated by multiple osillatory scalar fields and have shown that it 
works well at least in the case in which the scalar fields have incommensurable masses. In 
this formulation, we introduce an expansion parameter e that represents the ratio of the 
cosmic expansion rate to the masses of the scalar fields, and decompose the Hamiltonian 
is into an unperturbed part that depends only on the action variables and a perturbative 
part that is of the order e and bounded by a constant multiple of 1/t, where t is a time 
parameter of the system given by a 3 / 2 in terms of the cosmic scale factor a. Next, we look 
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for a canonical transformation that transforms the perturbative part to a quantity of higher 
order with respect to e and 1/t. As was shown in Ref. [8 ], we can construct such a canonical 
transformation iteratively, and by repeated applications of such transformations, we can 
transform the perturbative part to a quantity of the order of an arbitrary power of e and 
1/t. Furthermore, the new Hamiltonian system obtained by these transformations becomes 
solvable if the small perturbative part is neglected. By this method, we have proved that 
the Bardeen parameter is conserved with a good accuracy under the assumption that the 
masses of the scalar fields do not satisfy resonant relations. 

In the reheating phase, a dynamical instability caused by the parametric resonance 
plays a crucial role in the energy transfer from a macroscopic homogeneous mode to finite 
wavelength modes 9 ]. This instability can have a significant effect on the evolution of cos- 
mological perturbations [TTH 1131 ] . In fact, a numerical example showing a non-conservation 
of the Bardeen parameter was presented in ^] ] . In order to treat this problem, it is 
necessary to extend our Hamiltonian formulation to the resonant case. From this point of 
view, in this paper, we undertake this extension and with the help of it, we investigate the 
dynamical behavior of the universe dominated by multiple oscillatory scalar fields whose 
masses satisfy a resonance condition at least approximately. 

The present paper is organized as follows. First, in the next section, on the basis of the 
paper jU], we explain how to construct a solution to the perturbation equation in the long 
wavelength limit from an exactly spatially homogeneous perturbation of the background 
universe model. Then, in §3, we put the spatially homogeneous system of multiple scalar 
fields in an expanding universe into the Hamiltonian form and introduce the action-angle 
variables. 

In §4, for the case in which the masses of scalar fields satisfy resonance conditions 
approximately, we decompose the action-angle variables to fast changing variables and 
slowly changing variables, and reduce the dynamics of the system to that of slowly changing 
variables by a canonical transformation. Then, with the help of this formulation, we 
estimate the growth rate of perturbations of the system in the case in which the time 
parameter t of the system is smaller than 1/e. On the basis of this estimate and the 
analysis of some soluble examples, we argue that the Bardeen parameter of the system can 
grow if the system has a hyperbolic fixed point. Next, in §5, we analyze the evolution of a 
perturbation of the system in the time range t > 1/e and show that the Bardeen parameter 
is conserved with a good accuracy in this time range. §6 is devoted to discussions. In 
order to make the presentation clear, the proofs of most of the mathematical statements 
and the calculation of the growth rates of perturbations in concrete models are given the 
Appendices. 

Throughout the paper, the natural units c = h = 1 are adopted, and 8nG is denoted 
as k 2 . Further, the notation adopted in the articleP] is used for perturbation variables, 
and their definitions are sometimes omitted except for those newly defined in this paper. 
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§2 Evolution of cosmological perturbations in the long- 
wavelength limit 

As mentioned in the introduction, we can determine the dynamical behavior of cosmological 
perturbations on superhorizon scales by studying an exactly spatially homogeneous system. 
Since the perturbative analysis of the latter system can be used to determine its Lyapunov 
exponent, which is an index for a dynamical instability including chaos of a Hamiltonian 
system, we can also analyse the dynamical instability and integrability of the inflaton 
dynamics by such a study. 

In this section, we summarize the main results of the paper |U] in the case in which the 
universe is dominated by multi-component scalar fields, and explain how to construct a 
solution to the perturbation equation in the long-wavelength limit from an exactly homo- 
geneous perturbation of the model. We assume that the universe is spatially fiat (if = 0) 
throughout the paper. Hence, the background metric is given by 

ds 2 = -dr 2 + a(T) 2 dx 2 . (2.1) 

We consider the universe dominated by multi-component scalar fields whose energy-momentum 
tensor is given by 

Ty = Vy ■ V„0 - (V A • V A + 2U) . (2.2) 

In the long wavelength limit, the gauge-invariant variable Y{ representing the fluctuation 
of the scalar field 0; in the flat time slice can be expressed as 

Yi = 

W(X h X 2 ) := Xi ■ X2 — Xi ■ X 2 , (2.3) 

where \i is the combination of the exactly homogeneous perturbation of (pi, 5(pi, and the 
perturbation of the cosmic scale factor a, 5a, given by 

= Hi - 77 — • (2.4) 
H a 

From the equations of motion, it follows that this quantity always satisfies the equation 

nr^(l' x ) =const - (2 - 5) 

In general, the general solution to the homogenous dynamical system can be expressed 
in terms of the scale factor a and a set of integration constants as 

& = &(a,C). (2.6) 

For this expression, \i can De simpliy written as 

X. = (27) 
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As explained in the introduciton, one of the most important quantities describing the 
evolution behavior of perturbations is the Bardeen parameter defined by 



C = n - —a g , (2.8) 

where TZ and a g are perturbations of the three curvature and the shear of each constant 
time slice, respectively. In the present case, the Bardeen parameter ( can be written in 
terms of the gauge-invariant variable as 

6 ■ Y 

C = -H^—. (2.9) 

(0) 2 

Hence, from ()2.3|) . £ is represented in the long- wavelength limit as 

From (|2.5jl . the second term on the right-hand side of this equation is proportional to 
J dr/a 3 , which rapidly approaches a constant as the cosmic time r increases. Hence, in 
order to see whether the Bardeen parameter is conserved or not in the superhorizon stage, 
we can concentrate on the first term. For this reason, from this point, we assume that £ is 
expressed as 

f=-H^. (2.11) 

(0) 2 

§3 Evolution equations of corresponding exactly ho- 
mogeneous universe 

If we use the cosmic scale factor a as the time variable, the action for the homogeneous 
dynamical system introduced in the previous section can be put into the following Hamil- 
tonian form: 



S = / y^P^#i ~ hgda, (3.1) 
J i 

»° - T^l^EA + ^f' (3.2) 

i 

U((j>) = lY.rrtifi + Utett), (3-3) 

i 

where U- m t is assumed to be a sum of monomials in of degrees not less than 3. After 
changing the time variable to 
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this action can be expressed in terms of the non-dimensional canonical variables 



* = X"> ^HT^TTy (3-5) 



and the dimensionless parameters 



Hi = — , e = — K0o, (3.6) 
m 2 



as 



S' 



naQin (f)Q 



Now, let us introduce the action- angle variables (Ji,6i) by 



1/2 

f/ int (0o$) I d*(3.7) 



/2 1 
— VJi-cos^, (3i 



P^ = -V2^V^sin^, (3.9) 
which corresponds to the canonical transformation generated by the generating function 

W($,9,t) = -i^/i i t 2 $ J 2 tan^. (3.10) 

i 

Then, the Hamiltonian of the system is transformed to 

H = -C£^i + -fri^int) 172 - ^^sin2ft t . (3.11) 
Further, in terms of these new variable, the Bardeen parameter is expressed as 

C = ?7(^+4^nt) 1/2 ' 



3£ v ^ m\<g ' Ei/^i^(! -cos 20*) 
^(-<LPsin20* - P(l- cos 2e i )S6 i y (3.12) 



Evolution for t < lie 



In this section, we evaluate the growth rate of a perturbation of the system defined by 
the Hamiltonian (J3.ll)) during the period 1 < t < 1/e, assuming that the masses of scalar 
fields satisfy a resonant condition approximately. Before going to general arguments, we 
first explain the basic ideas by two simple models consisting of two scalar fields. 
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The first model is defined by 

£Ant(0) = Wife, 2/1x^2. (4.1) 

If we change the canonical variables from (6>j, Jj) (i = 0, 1) to (qi,Pi) (i = 0, 1) by the linear 
symplectic transformation 

0i = <?o #2 = 2g + <?i 

Ji = po - 2pi J 2 = pi, (4.2) 



we obtain 
where 

i/int can be written 



A*i^i + At2^2 = ^oPo + wiPi, (4.3) 
w = Aii ^i = -2/xi + /i 2 - (4.4) 



C/int = 7JiJ 2 1/2 {cos(2^-^) + cos(2^ + ^) + 2cos^} 



m^02 - — t 

= j(p - 2p!)pJ /2 {cosgi + cos (4g + q ± ) + 2cos(2g + ft)}, (4.5) 

where 

i = ^?-^- («) 



Hence, the contribution of this interaction term to the Hamiltonian is of order 1/t, and 
when averaged over q and q±, it becomes of order 1/t 2 . It generally holds on the interaction 
terms of the third degree of <fi. This feature will play an important role in the argument in 
§5. 

The second model is defined by 

tAnt(0) = \(/>l4>l, A*l ~ A*2- (4.7) 
By the linear symplectic transformation 

0\ = q 02 = qo + qi 

Ji=Po-Pi h= Pi, (4.8) 

we obtain 

V\J\ + H2J2 = u po + uip!, (4.9) 

where 

uj = pi uji = -Hi + (4-10) 
The interaction term is now written 

;U- mt = ^JiJ 2 {l + Jcos(20 1 -20 2 ) + ^cos(20i + 20 2 )+cos20 1 + cos20 2 } 



m§0§ t 2 1 z 1 ' 2 v 1 z > ' 2 

= S(Po-pi)pi 

X x 

{ 1 + - cos 2qi + - cos (4g + 2q-y) + cos 2q + cos (2g + 2ft) } , (4.11) 
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where 

/ii/i 2 

Hence, in this model, the contribution of this interaction term to the Hamiltonian is of 
order 1/t 2 , and its dominant part does not vanish even after it is averaged with respect to 
qo and q±, in contrast to the first model. 

In both of these models, lo\ represents the deviation from the exact resonance, which 
is given by 2/ii = /x 2 for the first model and by fi\ = yU 2 for the second model, respectively. 
Hence, lo\ is much smaller than u , and as a consequence, the motion of q is much faster 
than that of g x in general. For this reason, we call q and q± the fast angle variable and the 
slow angle variable, respectively. 

This behavior of the variables qi suggests that the dynamics of the slow variable is well 
described by a Hamiltonian H obtained from H by averaging it with respect to the fast 
angle variable q Q : 

1 f 2w 

H = — dq H. (4.13) 
2n Jo 

For the first model, this averaged Hamiltonian can be expressed as 

H=\{»- V) XI2 + \ \ J^q-M - ^)vT cos gi + 0(^1), (4.14) 
and for the second model, as 

H = -( U . p) 1 ' 2 + 31 1 (po - pi)pi{l + \ cos 2 Ql } + 0(^1), (4.15) 

where u> ■ p = u p + uipi . 

This simple procedure, however, does not give a useful approximation, since the higher- 
order terms with respect to 1/t still contains the large parameter 1/e, as is seen from the 
above expressions. 

Now, we will show for a generic system consisting of n (> 2) oscillatory scalar fields 
that this difficulty can be resolved by taking the average after applying an appropriate 
canonical transformation to the Hamiltonian. For that purpose, we transform the original 
action-angle variables to the fast canonical variables (<7o,Po) and the slow variables (<?i,Pi), 
as in the above examples. 

For the time being, suppose that the masses of the system satisfy n\ resonance relations 
of the form k fi = 0, where k is a vector with irreducible integer coefficients. Let R = (Rij) 
be a unimodular integral matrix such that its last n\ rows are given by the n\ vectors k 
defining the resonant relations, and consider the canonical transformation generated by 
W = PiRijOj, 

dW 

Qi = -^ = Rjk9 k , (4.16) 
dW 

Here, note that due to the unimodularity of R, we can assume that each of the new angle 
variables q^ also has the period 2n. 
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Now, let us decompose these new variables (q, p) into two sets of variables as 

(Qo)j = Qj, (Po)j =Pj (j = 1, • • • , no), (4.18) 

(<?l)i = Qno+i, (Pl)i = Pno+i (* = 1, ■ • • ) ( 4 - 19 ) 

where n + ri\ = n. Then, the set (qo,Po) becomes the fast variable and the set (<Zi,pi) 
becomes the slow variable, because p ■ J can be written as 

p-J = p -uj , (4.20) 

where 

(w )< := (* = " " ' > n o)- ( 4 - 21 ) 

In the above, we have imposed exact resonant conditions on the mass parameters. 
However, such strong conditions are rarely satisfied. In fact, it is known that for d > n — 1, 
the set of points in the u) space satisfying the condition 

inf \k\ d \(k ■ u)\ = (4.22) 

has measure zero, where 

\k\ := | fei| + |fe 2 | H h |fe„|. (4.23) 

Therefore, from this point, we only require that the mass parameters satisfy resonance 
conditions approximately. To be precise, we assume that 

inf \k\ d \k -w + k 1 -u 1 \=C>0, (4.24) 
fe ^o,feez n 

inf \k\ d \k ■ + fei ■ < eC, (4.25) 

where C is a constant of order unity. When u? satisfies f!4.24|) for some positive constants 
d, C, we say that u> is of the class Dq(c?, C). 

Next, we show that there exists a canonical transformation such that in a system 
obtained by that transformation, the dynamics of the slow variables can be determined 
independent of the fast angle variables with a good accuracy. For that purpose, first note 
that the Hamiltonian of the system can be written in terms of the variables qo,Po, Qi and 
Pi as 

H = -(u> ■ pa + u>i • Pi) 1/2 + A{q u p, t) + B(q,p, t), (4.26) 
e 

by Taylor expanding H with respect to 1/t, where 

(u>i); := Rnv+ijfij (i = 1, • • • , no). (4.27) 

We say that this Hamiltonian is of the type C m (a, Mi, M 2 , p), if the following conditions 
are satisfied for some positive constants a, Mi, M 2 , and p: 

(i) u> ■ p = ujq • Pq + u>i ■ p\ is bounded as 

\u-p\>a. (4.28) 
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(ii) tA can be extended to an analytic function in the domain -Di(p) in C n+ni+1 denned 
by 

D (p) ■= I (p,q, 7 ) ; Re P + l Im pl < p> l Im< ?l < p. 



Im- 



< Pj , (4-29) 
£>i(p):=£>(p)I»=o, (4.30) 



where D + p denotes the p-neighbourhood of an interval D, and in this domain, 
satisfies the inequality 

\tA\ < Mi. (4.31) 
Further, tA is periodic with respect to qi and real if (<?i,p, 1/t) are real. 

(iii) t m B can be extended to an analytic function in the domain D(p) C C 2n+1 defined 
above and satisfies the inequality 

\t m B\ < e m_1 M 2 . (4.32) 
Further, t m B is periodic with respect to q, real if 1/t) is real, and satisfies 

— I = a (4.33) 

Under this notation, the following proposition holds. 

Proposition 4.1 Let m be some positive interger, and consider the Hamitonian 
written in terms of fast canonical variables (q^ , Pq) and slow canonical variables (q[ m ^ , Pi"^) 
as 

= ^(u, • + «i • pS m) ) 1/2 + 4»(<li m) ,P (m) ,f) + B m (? (m) ,P (m, ,t). (4-34) 

Suppose that u) is of the class D (d,C) and that this Hamiltonian is of the type 
Cm(<J m , M[ m \ M^, p m ). Then, for any 5 > 0, there exists e > such that, for an 
arbitrary e satisfying 

< e < e , (4.35) 
there exists a function S m (q,p,t) satisfying the following conditions: 

(i) S m is periodic with respect to q and real if (q, p, t) is real. 

(ii) t m S m can be extended to an analytic function in the mult i- dimensional complex 
domain D(p m+ i), where p m+ i = p m — 5, and in this domain, satisfies the inequality 

(m) 



t m S m \ < e m LP (4.36) 



for some positive constant L\ 



M 
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(iii) Let H^ m+1 \q^ m+1 \ p( m+1 ), t) be a Hamitonian obtained from if( m ) by the canonical 
transformation generated by S m (q( m \p( m+1 \t): 

p (m) = p (m + l) + |^ (4 _ 3?) 



* {m+1) = * M + (4-38) 

= + ^ (4>39) 



= H( W . p (m+i))i/a + A m+1 ( g i m+i) ,p( m+i ),t) + B m+1 { q {m+1 \p {m+1 \mm 



e 

Then, if ( m+1 ) is of the type C m+ i(cr m+ i, M[ m+1 \ M^ m+1 \ p m +i) for some positive con- 
stants <Jm+i, M[ m+1 \ and M^"^ 1 " 1 , and the change of the A-term in the Hamiltonians 
satisfies the inequality 

\t m+1 {A m+1 ( qi ,p,t) -A m ( qi ,p,t)}\ < yM< m+1) , (4.41) 

for {q!,p, l/t) E Di(p m+ i). 

(For the proof, see the appendix A). 

Since the generating function S m in this proposition is of the order e m /t m , we have B m ~ 
e m_1 /t m and B m+ i ~ e m /t m+1 . This implies that we can make the part of the Hamiltonian 
that depends on the fast variables q$ aribitrarily small by taking the original Hamiltonian as 
the starting point (m = 1) and applying canonical tranformations given in the proposition 
repeatedly. Therefore, we can expect that the evolution of the slow variables (<fi,Pi) can 
be determined with a quite good accuracy from an effective Hamiltonian H( m \q[ m \ p^ m \ t) 
obtained by dicarding B m term from H {m \ if we take m sufficiently large. In this trancated 
system, p^ becomes constant, and the behavior of q^ can be obtained by a simple time 
integration of a function that has a definite t-dependence when p^ and a corresponding 
solution for (q[ , p[ m ') are given. The evolution of the original variables can be calculated 
from this solution by applying the known canonical transformation connecting these two 
sets of variables. 

Next, we show that this expectation is correct by estimating the errors produced by the 
trunction. Let us use the symbol A to represent the difference of a quantity for the exact 
system and the corresponding quantity for the truncated system. Further, let us denote 
a quantity for the truncated system by the upper case letter of the lower case symbol 
representing the corresponding quantity for the exact system. For example, the errors of 
canonical variables are denoted as 

AQ = q - Q, (4.42) 
AP = p - P, (4.43) 

and the errors of solutions to the perturbation equations for both systems are denoted as 

A5Q = 5q- 5Q, (4.44) 
A5P = 5p- 5P. (4.45) 
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Further, let us denote the sets of slow variables, (qi,Pi) and (Qi, Pi), by z and Z, respec- 
tively. Finally, for a function f(t), let us define ||/||(t) by 

||/||(f) = sup |/( S )|. (4.46) 

l<s<t 

Before evaluating AZ, A5Z and SZ, we give the general technique to evaluate an upper 
bound on the norm of the solution X to the first order differential equation: 

jX = VIX + S, (4.47) 

where X, S are N column vectors and Q is N x N matrix. For the detail, see the appendix 
of the paper [3]. If we define the norm of the solution X by 

\X\ 2 :=X j X, (4.48) 

we decompose Q into a sum of Q%, Q2' ^i from the perturbed part is of order 1/t and S^2 
from the unperturbed part is of order e and A m j is the largest eigenvalue of the hermitian 
matrix Qh% '■— + \X\ is bounded as 

| W a m+ | S |, (4.49) 

where 

A := Ami + A m2 . (4.50) 

If A is an eigenvalue of Sim, —A is also an eigenvalue of £lm, because fij can be written in 
the form of I Si where Si is the symmetric matrix and I is defined by f)4.79j) . So A m j are 
non-negative. Therefore the norm of the solution is bounded as 

A 



\X\ < exp 



■dt 

2 



|X(1)| + / dt\S 



(4.51) 



ui 

In evaluation of the contribution from the source term, we use 

k k 

I <Z)l«il» ( 4 - 52 ) 

i=i i=i 

which is obtained from the Cauchy Schwarz inequality 

\A ] B\ < \A\\B\. (4.53) 
Since the different definition of the norm 

\X\ m := max \(X)i\ (4.54) 

l<i<N 

satisfies 

1 \X\ < \X\ m < \X\, (4.55) 



we can identify \X\ and \X\ m assuming that the number of degrees of freedom is of 
order unity, so we simply omit the subscripts m. 

Then, the truncation error for the fourth-order system can be estimated as follows. 
From now on we omit constant coefficients of order unity except V. 



11 



Proposition A.2A Let T be an upper bound of the eigenvalues of the the hermitian 
matrix tflni/^ defined by 



/ d 2 A 



fii : = 



d 2 A \ 



dQ 1 dP 1 

d 2 A 



dP 1 dP 1 

d 2 A 



(4.56) 



\ dQ 1 dQ 1 dP 1 dQ 1 / 



where A is the A-term of the fourth-order Hamiltonian H = and all the elements of 
Qi are bounded by /5 2 t in the domain Di(p^ — S). This quantity T gives an upper 
bound on the growth rates of the errors AQi, APi, ASQi, and AS Pi, and the growth of 
the perturbations SQi, 5 Pi. 
Let us define (3 by 

P 



1 



r + i' 

Then, for an arbitrary positive e, in the time interval 

1 



1< t < 



the truncation errors of canonical variables are given by 

|AQo 



|AP | < e'\ \AZ\ < e% 

and the errors of perturbation variables are 

ASQq 



t 



< 



(4.57) 



(4.58) 



(4.59) 



< e 2 |5Q (l)| + e||(5Po||+e||5Z| 



\\ASP \\ < e 3 \SQ (l)\ + e 2 \\SP \\+e ;i \\SZ\ 
\\ASZ\\ < e 2 \SQ (l)\ + e\\SP \\ +e\\SZ\\, 

under the initial conditions: 



(4.60) 



APo(l) = AZ(1) = AQo(l) = 0, (4.61) 
ASQ (1) = A5P (1) = A5Z(1) = 0. (4.62) 

(For the proof see Appendix B.l.) 

This proposition shows that the truncation errors can be made small in the fourth-order 
system (m — 4). In order to obtain the information on the behavior of the original variables, 
we also have to estimate how this truncation error affects the original variables through the 
canonical transformation which connects the fourth-order system and the original system. 
The next proposition gives that estimate. 
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g(4) )p (4)) by the 



Proposition 4.25 The difference between (q^,pW) obtained from 
canonical transformation and (Q^\ P^) obtained from (Q^\ P^) by the same canonical 
transformation has the upper bound 



AP n (1) |<e 3 , |AZ (1) | < 



(4.63) 



and the corresponding difference in the perturbation variables is estimated as 



A5Q 



(i) 



< e 2 \6Q$\l)\+e\\6pW\\ + e\\6ZW\\, 



(4), 



A5P (1) || < e 2 |5Q^(l)|+e||5P c 



(4)/ 



\ASZ^\\ < e 2 |5Qo(l)| +e||5P c 



>(4)| 
>(4)| 



'\\5Z^\\. 
\5Z^\\, 



(4.64) 



in the same time interval as in the previous proposition. (For the proof see Appendix B.2.) 

These two propositions show that the truncation errors can be made small even with 
respect to the original variables, if we truncate the system at the fourth-order. This order 
is minimal in the sense that the trunction at a lower-order system produces errors of order 
unity in the perturbation variables. Conversely, if we go to higher-order systems, we can 
obtain a better approximation. Further, we can prolong the time interval in which the 
approximation is good. In fact, we can show that in the m-th order system with m > 4, 
the same estimates for the truncation errors as in the above propositions hold in the interval 



1< t < 



1 



,(m-3)/(r+l) 



(4.65) 



In particular, for m larger than T+4, the approximation is good in the interval 1 < t < 1/e. 

Finally, let us evaluate the growth of perturbations of a truncated higher-order system. 
The Hamiltonian of such a system can be in general expressed as 



where 



H = -(w • P + wi • Pi) 1/2 + A(Q 1 , P, t), 



\tA\ < M, u>! = 0(e). 



(4.66) 



(4.67) 



As mentioned before, the fast action variables Pq for this system become constants of 
motion, and the equations of motion for the slow variables (Pi, Qi) do not contain small 
parameter e essentially. Therefore, in contrast to Po, the slow action variables Pi are not 
conserved in general, and as a consequence, the Bardeen parameter can grow considerably 
in the time interval 1 < t < 1/e. (For estimate of the growth rates of the Bardeen parameter 
in the concrete examples, see the Appendix E.) This can be confirmed by the following 
estimate of the upper bound of the perturbation variables 5Q and 5P. 
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Proposition 4.3 In the time interval 1 < t < 1/e, the perturbation variables SQ and 5P 
are bounded as 



\5Z\ < t r [\5Z(l)\ + \5P (l)\(t-l)], 
\5P \ = \SP (1)\, 

\SQ \ < \8Q (l)\ + -J5P (l)\(t-l) 

t r+1 - 1 f +2 - 1 t r+1 - 1 



r + i 



r + 2 



r + i 



(4.68) 
(4.69) 



(4.70) 



where T is the quantity defined in the proposition 4.2B. In particular, the growth rate of 
the perturbation variables is not exponential and at most a power of t. 

proof The evolution of the slow variables SQi and 5 Pi is determined by the equation 

/ d 2 A 



d_ f SQi 
dt V SPi 



uji ■ sjJi d 2 A \ 

dQidPi ~2e(u -P) 3 / 2 + dPidPi 
d 2 A d 2 A 



V dQidQi 
+SP 



+ 



dPidQi 

d 2 A \ 



SQi 
6Pi 



2e(w • Pf/ 2 dP Q dPi 
d 2 A 



\ 



dP dQi / 
By evaluating the coefficients in the right-hand side, we obtain 

±\5Z\<(j + e)\5Z\ + \5P (l)\, 

where we have used 5P = 5P (1). By integrating this inequality, we obtain 

\5Z\ < t T exp[e(t - 1)}[\SZ{1)\ + \5P (l)\{t - 1)]. 
Since we are considering the time interval 1 < t < 1/e, we have 

\5Z\ <t r [\5Z(l)\ + \5P (l)\(t-l)]. 
Next, 5Q obeys the equation 

1 



dt ^° ~ 2 e (w • Pf' 2 



(w • 5P + wi • SPi) 



, tn d ^ d .„ 9 , dA 



dP ' ^ dQi ' 1 dPi'dPo' 
Integrating the right-hand side of this equation with respect to t yields 



(4.71) 



(4.72) 



(4.73) 



(4.74) 



(4.75) 



\SQ \ < \5Qo(l)\ + - e ^dt\5P (l)\ + ^dt\5Pi\ + J^[\6P (l)\ + \6Z\] 

< \5Q (1)\ + ^fdt\5P (l)\ + J^dt\5Z\. (4.76) 
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From (|4.74|) and this equation, we obtain ()4.70|) . □ 

Next, we show that the above general estimate of the upper bound on the growth of 
perturbations is rather good, by presenting an example in which the upper bound is nearly 
saturated. Let us consider a system in which the Hamiltonian flow in the phase space 
(Qi,Pi) has a equilibrium fixed point. For simplicity, we assume the exact resonance 
condition u>i = 0. Then, in terms of the phase space variable Z defined by 

Z i = Q\ (1 < % < m), Z ni+i = P{ (l<i< m), (4.77) 

the Hamiltonian equations of motion can be written 

dZ 8H 

M = ' flZ" < 478 > 

where / is the matrix of degree 2n\ expressed in terms of the unit matrix E of degree n\ 
as 

' £ N 



We consider the Hamiltonian of the form 

H = ■ P ) 1/2 + ^a(Q u P), (4.80) 

where 7 = 1 or 7 = 2. When the resonant interactions come from interaction terms of 
the third order in the scalar fields 0, 7 = 1, and when they come from interaction terms 
of the fourth order in 0, 7 = 2. We focus on the dynamical behavior near the fixed point 
(Po(l), Z{1)) and expand a(P , Z) as 

a = a + b T Po + \P^CP, + Z t DPq + \z T FZ, (4.81) 

dropping the terms of degree not less than 3 with respect to deviations from the equilibrium; 
Pq = Pq — Po(l), Z_ = Z — Z(l), where T implies to take the transposition of a matrix, 
ao is a constant, b and Pq are no-dimensional vectors, C is an n x n symmetric matrix, 
D is a 2ni x no matrix, and F is a 2n\ x 2ni symmetric matrix. Since Z(l) is the fixed 
points, a does not contain g T Z^. Then, the evolution equations of the slow variables Z can 
be written in the matrix form as 

f = L[IFZ + IDBi. (4-82) 
Taking the variation of this equation and using the relations 

oP = oPo(l), (4.83) 

we obtain 

j t 5Z = ^[IF6Z + ID5P Q (1)]. (4.84) 
Here, note that the 2ni x 2ni matrix X := IF satisfies the equation 

X T = IX I. (4.85) 
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From this and I 2 = 



— 1, it follows that the characteristic polynomial 



A(A) = det(X - \E) (4.86) 

is an even function of A: 

A(A) = A(-A). (4.87) 
In addition, since A(A) is real polynomial, we have 

A(A) = A(A*). (4.88) 

Hence, if a + bi, (a, b real) is an eigenvalue of X, all of ±a ± bi are also eigenvalues of 
X. From this, it follows that if X has an eigenvalue with a non- vanishing real part, the 
fixed point is unstable. We say that such a fixed point is hyperbolic. In constrast, if all 
eigenvalues of X are pure imaginary, the fixed point is stable. We say that such a fixed 
point is elliptic. In the case in which the fixed point is neither hyperbolic nor elliptic, 
i.e., some of the eigenvalues of X vanish and the other eigenvalues are pure imaginary, 
the flow around the fixed point is stable in the linear analysis but may become unstable if 
higher-order terms are taken into account. 

To be precise, this stability argument applies to this autonomous systems. In the present 
case, due to the existence of the time- dependent factor 1/t 7 , the real stability depends on 
the value of 7. Since the system under consideration is linear, we can check it directly by 
solving the equation. First, by diagonalizing the matrix X as 

S^XS = A, (4.89) 

where A is a diagonal matrix, the general solution for SZ can be written 

SZ = Sexp[ / -AjS^K - X- l ID5P Q {\), (4.90) 
Ji t 1 

where 

K = 5Z(1)+X- 1 ID5P (1). (4.91) 
Next, the equation for SQo is given by 

Hence, SQ can be expressed as 

,Q = ,Q 0(1) -g -^W- ( f - 1)+B , (4.93) 

where 

R = C5P (l)j^ + l^5Z T D. (4.94) 
When 7 = 1, inserting the above expression for SZ into this equation, we obtain 

t x - 1 

|#| < — — + lnt, (4.95) 
A 
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for A > 0, and 



(4.96) 



for A = 0, where A is the maximum of the real parts of the eigenvalues of X. Thus, 5Qo is 
always unstable. In contrast, when 7 = 2, \R\ is bounded from above as 

\R\ < e A + 1, (4.97) 

and the stability of 5Qo depends on the value of 5Po. 

Finally, we give an exactly soluble example such that there appear both hyperbolic 
fixed points and elliptic fixed points and the phase flow pattern is similar to that of the 
pendulum in a conservative field. It is given by the Hamiltonian 

H = ^> Po) 1/2 + ^ QaP* + fccosQxJ + ^, (4.98) 

where a, b, c are positive functions of P . 
Since there exists a conserved quantity, 

C = + focosQi = const, (4.99) 

in this model, we can determine the phase flow pattern in the (Qi,Pi) plane easily. In 
particular, we find that the phase plain is divided into a region of oscillatory motions and 
a region of rotationary motions by the separatrix defined by 




(4.100) 



Further, we find two types of fixed points. One is elliptic points, Pi = 0, Qi — (2k + 1)tt, 
where k is an integer, each of which is surrounded by flows corresponding to oscillatory 
motions. The other is hyperbolic fixed points, Pi — 0, Q\ — 27ck, where k is an integer. 
These fixed points are connected by the heteroclinic orbits (|4.100j) . This implies that near 
the heteroclinic orbits, the perturbations 5Q and 5P, therefore the Bardeen parameter 
grows. 



§5 Evolution for t > 1/e 

We have treated the case in which uj is of the class Do(d,C) and satisfies the resonance 
relations at least approximately (|4.25l ) for t < 1/e. When there exists deviation from the 
exact resonance, as the time proceeds, the instability characteristic of resonance disappears 
and the evolution of the system is reduced to the superposition of oscillations with ampli- 
tudes and frequencies of different orders. We demonstrate this fact by investigating the 
case in which u> is of the classes Do(d, C) and D(d, eC) for t > 1/e: we say that u) is of the 
class D(d, C) if 

inf \k\ d \k ■ uj + k t ■ ujxl = C, (5.1) 
is satisfied for some positive constants d, C. 
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By redefining t, H by et, H/e respectively, we investigate the evolutionary behavior for 
t > 1 of the system defined by the Hamiltonian of the type Ci^(a, M n , M 12 , M 2 , p): we 
say that this Hamiltonian 



2 



H = ^(uJ-p) 1 / 2 + A(q 1 ,p,t) + B(q,p,t), (5.2) 



e 



A(g ljP ,t) = D(p,t) + E(q u p,t), (5.3) 

is of the type C ni: i(a, M u , M 12 , M 2 , p), if the following conditions are satisfied for some 
positive constants a, M n , M 12 , M 2 , and p: 

(i) cl> ■ p = ujq ■ p + u?i ■ pi is bounded as 

|u> • p| > a. (5.4) 

(ii) t 2 D can be extended to an analytic function in the domain D 2 (p) := D(p)\ q=0 in 
C n+1 and satisfies the inequality 

\t 2 D\ < eM n . (5.5) 
Further, t 2 D is real if (p, l/t) are real. 

(iii) t'E can be extended to an analytic function in the domain Di(p) C C n+ni+1 and 
satisfies the inequality 

\t l E\ < e' _1 M 12 . (5.6) 
Further, t'i? is periodic with respect to q 1; real if (<7i,p, l/t) is real, and satisfies 

27T /-27T 



d ni q x t l E = 0. (5.7) 



(iv) t m £? can be extended to an analytic function in the domain D(p) C C 2ri+1 and satisfies 
the inequality 

\t m B\ < e 2(m - 1} M 2 . (5.8) 
Further, t m _B is periodic with respect to q, real if (q, p, l/t) is real, and satisfies 

-27T ^271 



(27T) 



/ • • • / d no q t m B = 0. (5.9) 
Jo Jo 



Since the interaction terms of the third degree of do not contribute D, D begins with 
terms of order e/t 2 . 

First we show that we can make q$ depenent part B small by the canonical transfor- 
mation S. We say that the Hamiltonian is of the type C^(<7, M 1: M 2) p) if the conditions 
with e replaced with e 2 in the definition of the type C m (a, M 1; M 2 , p) are satisfied. 
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Proposition 5.1 Let m be some positive interger, and consider the Hamitonian 
written in terms of fast canonical variables (q^ , Pq) and slow canonical variables (q^ , f>["^) 

as 

# (m) = |(«o -pi n) + Wl - P s m) ) i/2 + A m ( 9 ; m) , P ( m ),t) + 5 m ( 9 (™),p( m ),t). (5.10) 

Suppose that a; is of the class Do(d,C) and that this Hamiltonian is of the type 
C^(a m , M[ m \ p m ). Then, for any 5 > 0, there exists e > such that, for an 

arbitrary e satisfying 

0<e<e , (5.11) 
there exists a function S m (q,p,t) satisfying the following conditions: 

(i) S m is periodic with respect to q and real if (q, p, t) is real. 

(ii) t m S m can be extended to an analytic function in the multi-dimensional complex 
domain D(p m+ i), where p m+ i = p m — 5, and in this domain, satisfies the inequality 



t m S m \ < e 2m Li m) (5.12) 



for some positive constant L^™\ 



(iii) Let H tym+l \q^ n+1 \p^ m+l \ t) be a Hamitonian obtained from by the canonical 

transformation generated by S m (q( m \ p( m+1 \ t): 

P M =P (m+1) + ^, (5-13) 

1 l '" +t > = 8 <m> + g^I) • 

ff<"- + " = H lm) + ^ (5.15) 
at 

= - 2 (u,.p^Y 2 + ^ m+ i(qi m+1) ,P (m+1) ,t) + B m+l {q {m+1 \p {m+1) M^) 

Then, if ( m+1 ) is of the type C% l+1 (a m+ i, M[ m+l \ M^ 71 * 1 ^, p m+ i) for some positive con- 
stants a m+ i, M[ m+1 \ and M^™ 4 ^, and the change of the A-term in the Hamiltonians 
satisfies the inequality 

2m 

\t m+1 {A m+1 ( qi ,p,t) - A m (qi,P,t)}\ < — M 2 (m+1) , (5.17) 

for (q u p, 1/t) G L>i(p m +i). 

The proof is obtained by replacing e with e 2 in the proof of the Proposition 4.1. 

We show that the fast action variables po are perpetually stable and oscillate around 
the initial values with amplitudes of order e 2 and with frequencies of order 1/e 2 . 

Proposition 5.2 There exists some constant C such that 

IpP - p£\l)\ < e 2 C . (5.18) 
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proof We consider the transformed Hamiltonian given by 



H^ = -J^.p^fl 2 + A 2 + B 2 , (5.19) 



where 



\tA 2 \ < M[ 2) , (5.20) 
\t 2 B 2 \ < e 2 M 2 (2) (5.21) 



which gives 



We obtain 



dp£ ] 8B 2 



dt " d q {2) ' 

.2 ^(2) »A2) 



(5.22) 



By using the transformation law 

AC c 2 T-(l) 



we obtain 



\vf-vfiX)\ 

< - Pfl + Ip?' - '(i)l + Ip^ 2) (i) - 

'(1) M (2) r(l) 

J _L , 2 1V1 2 , 2 -^1 



< e 2^ + e 2^ + £ 2^ ) (5>25) 



which completes the proof. □ 

Next we show that we can make q\ dependent part E small by the canonical transfor- 
mation T. 

Proposition 5.3 Let m, / be some positive intergers, and consider the Hamitonian H^ m ' 1 ^ 
written in terms of fast canonical variables (<?(j m Po™ ''') an d s l° w canonical variables 
{q^Jr*) as 

# (m '° = |(u -p^Wpi m 'V 2 +^ 

(5.26) 

Suppose that uj is of the class D(d,eC) and that this Hamiltonian is of the type 
C m ,i{cm,i, My^' l \ M^ ,l \ M 2 m ' l \ p m j). Then, for any 5 > 0, there exists e > such that, 
for an arbitrary e satisfying 

< e < e , (5.27) 
there exists a function T m j(qi,p,t) satisfying the following conditions: 

(i) T m j is periodic with respect to qi and real if (qi,p,t) is real. 
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(ii) t l T m j can be extended to an analytic function in the mult i- dimensional complex 
domain D 1 (p m i+1 ), where p m ,i+i = p m ,i—S, and in this domain, satisfies the inequality 

\t l T m!l \ <e l L { 2 m ' l) (5.28) 

for some positive constant L 2 m,/ ' ) . 

(iii) Let H ( - m ' l+1 \q ( - m ' l+1 \p( m ' l+1 \t) be a Hamitonian obtained from H^ m,v> by the canon- 
ical transformation generated by T mi ;(qr| m, '' ) , p( m >' +1 ), £): 



p(^.0 =p W+i) > (5.29) 



n (m,l+l) _ (m,l) , dT m,l /, oi \ 

#(m,I+l) = #(m,0 + ggW (5.32) 

= |( W . p W+'))^ + D mil+1 ( P KW),() + S mi , +1 (q^ +1 ), p(^+D, *) 
+fl m , + i(g ( ^ ,+1 \p (m ' +1) ,f) (5.33) 

Then, H^ l +V is of the type C m , l+1 (a m , l+1 , m[?' 1+1 \ m[^ 1+1 \ M 2 (m ' m) , p m , m ) for 
some positive constants <7 m ,H-i> M-f™''" 1 " 1 ^, M-[™' i+1 \ and M^ m, ' +1 \ and the change of 
the D-term in the Hamiltonians satisfies the inequality 

\t l+1 {D m , l+1 (p,t) - D nhl (p,t)}\ < ^M^ l+1 \ (5.34) 

for (p, 1/f) e D 2 {p m ,i + i). 

{For the proof, see the appendix C). 

We show that the slow action variables pi are perpetually stable and oscillate around 
the initial values with amplitudes of order e and with frequencies of order 1/e. Notice that 
T is independent of qr . 

Proposition 5.4 There exists some constant C\ such that 

\pf -pf^^eC^. (5.35) 

proof We consider the transformed Hamiltonian H^ 171 ^ {m = 2, I = 2) given by 

tf(^) = 1 ( W . p P,2))i/2 + £, 22 + E22 + S22) (5. 36) 



where 



t 2 Ad 


< 




\t 2 E22\ 


< 


eMS 2 \ 


\t 2 B22\ 


< 


e 2 M 2 (2 ' 2) 



(5.37) 
(5.38) 
(5.39) 
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The evolution equations for the slow action variables pi are given by 

dpf' 2) dE 22 dB 22 e 



dt " dqp 2) dq^ 2) ^ 

which yields 



(5.40) 



\p {2 > 2) -pf 2 \l)\<eC, (5.41) 
for some positive constant C. 

\p^\t)-p^\l)\ 

< \p^\t)-p™(t)\ + \pM\t) -p?' 2 \l)\ + |p< 2 ' 2 >(l) -p^(l)| 

< eCi. (5.42) 

□ 

We investigate the dynamical behavior of the cosmological perturbations. We obtain 
the Hamiltonian of the type £3,3(03,3, M^' 3) , M$ 3 \ M 2 (3 ' 3) , p 3)3 ) by carrying out the 

canonical transformations in the order Si, S 2 , T 31 , T 3 2 . 

Later we omit the superscript (3,3) and constant coefficients of order 1. 

From now on we assume that 

|<Jqfo(l)| ~ l^o(l)| ~ \5qi(l)\ ~ \5 Pl (l)\ ~ 1, (5.43) 

since in the linear perturbation, the scale of the perturbation variables is arbitrary. Then 
we obtain the proposition below. 

Proposition 5.5 The transformed perturbation variables are estimated as 

\5p -5po(l)\ <e 2 , (5.44) 
|<fpi-<fpi(l)| <e, (5.45) 

Sqo = -~^(t - 1) . p | 1))3 / 2 (^ • Spo(l) + on • 5 Pl (l)) + i2 , (5.46) 

^ = "57^ ~ 1} (^-p(l))3/2 ^° ' + ^ ' + Hl ' ( 5 - 4? ) 

where the residual part Rq, Ri is bounded as 

\Ro\ < Co(t-l), (5.48) 

|-Ri| < eCi(*-l), (5.49) 

for some positive constants C , Ci. 

For the proof see the appendix (subsection D.l). 

We have evaluated the canonically transformed variables 5q^ 3 ' 3 \ op( 3,3 ). Then we eval- 
uate the difference between such transformed variables and the original variables Sq^ 1 ' 1 ^, 
<jp(i,i). 
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Proposition 5.6 



1^(1,1) _ Jg (3,3) | < 1; 

\Sp^)-6p^\ < 1. 



(5.50) 
(5.51) 



For the proof see the appendix (subsection D.2). 

According to propositions 5.5, 5.6, for t > l/e, the Bardeen parameter ( stays constant 
in a good accuracy. 

§6 Discussion 

In this paper we have constructed the method for analyzing the resonance phenomena 
of the Hamiltonian system obtained from the multiple oscillatory scalar fields in the ex- 
panding universe. We have shown that the perturbations including the Bardeen parameter 
can grow when the resonant interactions between the homogeneous modes exist. If the 
truncated Hamiltonian system has hyperbolic fixed points in the phase space of the slow 
variables, the perturbations grow at the speed of a power of t in the three leg interaction 
systems, while the growth of perturbations is bounded from above in no less than four leg 
interaction systems. In the models where u?x = 0, we have found the network constructed 
by the heteroclinic orbits around which perturbations are unstable. On this network, the 
orbits become irregular, unperiodic, complicated and probabilistic. In order to evaluate 
the amplitude of the Bardeen parameter which is directly related the cosmic structure 
formation, we must calculate the measure occupied by the stochastic network in the whole 
phase space of the homogeneous modes. 

Recentlty the relaxation in reheating is investigated numerically and is interpreted in 
terms of the turbulence phenomena of the dynamical system of large number of degrees of 
freedom [IH]. In order to understand the turbulence phenomena, our fixed point analysis 
of the Hamiltonian system may be useful. 
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§A Appendix; Proof of Proposition 4.1 

We consider the canonical transformation induced by the generating function given by 




(A.l) 



feo^o 
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where k ■ q( m "> = fc • q^ + fci ■ q[ m ^ and the sum is taken over k such that k ^ 0. B m is 
decomposed as 

B m (p {m+1 \q (m \t) = J2 b k (p( m+1 \ty k -* im \ (A.2) 

fco^o 

where the sum is taken over k satisfying k ^ 0. 
The transformed Hamiltonian is 

jj(m+l) _ jj(m) _|_ UkJ m 

dt 

_2, < m+ i )V/ 2 i 1 1 dS m (q {m \p {m+1) ,t) 
--^■P ) + ~ (w . pKi))!^ + 

+A m (g! m+1) ,P (m+1) ,t) + J R 2 

no f„(m) „(m+l) j.\ 

+ fi m ( q M, P ^),t) + R 3 + — (A.3) 

where ■ ■ ■ R3 are defined as 

* = -P (m) ) V2 " *(« "P (m+1) ) 1/2 " 7 (u? . p( l 1))1/2 ' (P (m) (A.4) 

i? 2 = A m (qi m) ,p("V) - A m (g! m+1) ,p( m+1 ),t), (A.5) 
i? 3 = B m (q^\p^\t) - B m (q {m \p (m+1 \t). (A.6) 

We determine the generating function S m so that the leading term depending on the fast 
angle variables q$ can be eliminated; 

1 1 . dS m (qWpWt) _ (m) (m+1) 

By comparing the Fourier components in the both hand sides, we obtain 

S k ( P (m+1 \t) = ic ( w .p("H-i))i/=»^_ 6fc (p(™+i) >t ). (A.8) 

(u> ■ k) 

Before we prove the analyticity and evaluate the upper bound of S m (q( m \ p( m+1 \ i), we 
use the following lemma. 

Lemma 

As for the Fourier series 

F(q) = J2 F ^ kq > ( A -9) 
fc 

where 

n 

k-q = J2 k i<H- ( A - 10 ) 

i=i 

(1) If F(q) is analytic and satisfies 

\F(q)\<C, (A.ll) 
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on the domain |Imq| < p, then 

\F k \<Ce-\ k \ p , (A.12) 

where 

|fe| = life! | H h \k n \. (A.13) 

(2) If on the domain |Imqr| < p, F k satrisfies 

\F k \ < Ce~ wp , (A.14) 

then F(q) is analytic on the domain |Imqr| < p, and for an arbitrary 5 (0 < 5 < p), on the 
domain |Imqr| < p — 5 

\F(q)\ < -p, (A.15) 

where we assume 5 < 3. 
Proof of Lemma 

(1) The Fourier coefficients F k are given by 

-i /*27r /*27T 

By Cauchy's theorem, the path of integration in the above integral can be shifted to 
qi = xi ± ip, < Xi < 2n where we choose the sign equal to — fcj, we get 

|F fc | < Ce~ wp . (A.17) 

(2) For an arbitrary positive 5, (0 < 5 < p), on the domain |Imq| < p — 5, 

\e ik ' q \ < e |fe|(p - 5) . (A.18) 



On the domain |Imq| < p — 5 



\F\ < ^2 \ F > 



k e 



ifc-q I 



< ^ Ce -I% e |fe|(^) 
fc 

k 

= C(l + 2^e- fc5 ) n 

fc; 

1 + e" 



fc>0 
-5 X " 



c 



< cA n , (A.19) 



where we use the inequality 



1 - e 

1 

S 

1 + e~ & i 
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< (t)> (A.20) 



for < 5 < 3. Namely on the domain |Imqr| < p — 8, F(q) converges uniformly and 
absolutely, therefore F(q) is analytic. Since 5 is arbitrary, F(q) is analytic on the domain 
|Imq| < p. 

Proof End of Lemma 

We evaluate the right hand side of (|A.8|) . For arbitrary p( m+1 ) e D(p m ), we put r m as 
the maximum value of u> ■ p( m+1 ) ; 

o m < w ■ p (m+1) < r m . (A.21) 
Since t m B m is analytic and bounded by e m_1 M2 on the domain D(p m ), 

\t m b k {p {m+l \t)\ < e m_1 M^ m) e- |fc|Pm . (A.22) 

As \k\ increases, while by, decay exponentially, the contribution from the denominator grows 
like power; 

1 \k\ d 

The exponential decay defeats the power grow. We use the inequality 

[jfel* < e W*(__)<*_L (A.24) 
e o d 

for an arbitrary positive 5. The proof is as follows. f(x) = x — dlnx has a minimum at 
x = d. Therefore 

31 S ? - (A.25) 



For x = |fe|5, we obtain this inequality. We evaluate the Fourier components Sk as 
| r5fe(p(m+ i) )f) | < e r r y 2 i^e m - 1 Mf ) exp(-|fc|p m ) 



C 

M (m) d 1 

< e m rV 2 ^(-) rf -exp{-|fc|(p m -5)}. (A.26) 



Therefore t m S" m is analytic on the domain D(p m ) and bounded as 

- (m) 



\t m S m \ <e m Lr\ (A.27) 



where 

4 m '=^^(^, (A.28) 

on the domain D(p m — 25) for an arbitrary positive 5 < 3. We put p m +i = Pm — 35 and on 
the domain D(p m+1 ) 

ao rim) 

| p M_ p (tn+i)| < |3^)|<e m ^(l + Pm+i r, (A.29) 

|„(m+l) (m)| ^ I i < m^l /-. , \m / A QfA 

19 ~<? I < I 9p (m+l) I - 6 -^-(i+Pm+lJ • (A.3Uj 
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Therefore if e satisfies 



e m ^_(l + Pm+1 ) m <35, 



when (p^+D^+D) e £>(p m+1 ), (p("0,gM) E D(p m ). 

Next on the domain D(p m+ i), we evaluate the residual part R defined as 



R — R\ + -R2 + -R3 + 
For evaluation, we use the Taylor expansion; 



ds m (q (m) ,p (m+1 \t) 



dt 



n-l 

l^(i)-^(o)-E^ F(fc) (°)l^^[ F(n) ]' 



k=l 

where [A] means the maximum value of \A\. As the function we take 

F(Z)=f(x + Z(x-x )), 
where x is the multidimensional vector. 

df 



F«(0) 

F (2)(Q) 

The residual part is 
As for i?!, 



dx i 



(x - x )\ 



d 2 f 



dx i dx^ 



(x - x y(x - x y. 



dx 11 ■ ■ ■ dx %n 



(x - xqY 1 ■■■(x- x y 



1 UJiUJj (p (m) - p( m+1 )) ! (pH - p( m +!))i 



4e (u> • p) 3 / 2 



Then we evaluate 



t 2m RA < 



1 1 J e m L\ 



(m) 



4e (x 3 / 2 



um- 



As for i? 2 , 



life I < 



Then we evaluate 



As for i? 3 , 



Ji/f M r M 

|r +1 i? 2 | < 2n^— e m ^-. 





|i? 3 | < 



dp 



(A.31) 



(A.32) 

(A.33) 

(A.34) 

(A.35) 
(A.36) 



(A.37) 

(A.38) 
(A.39) 

(A.40) 

(A.41) 
(A.42) 
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Then we evaluate 

\t 2m R 3 \ <ne m ~ x -^— e m =V. (A.43) 
o o 

In addition, 

|r +1 ^| < e m L { r\m + 1+pm+1 ). (A.44) 

Therefore we obtain 

\t m+1 R\ < e m M, (A.45) 

where 



M = ^(l+^r 1 1 



4^ 2 



A' 
um — : 



(m) 



2 



+ 2n- 1 



5 2 



ii j(m) 7- (m) , 

+6-^(1 + p^)- 1 ^ 2 ^ 1 + 4 m) (m + i±|^). (A.46) 

We decompose R into the slowly varying part R s and the fast varying part Rf as 

R = R s + R f , (A.47) 

where 

^ = ^i? ofel (p (m+1) ,Oexp[ife 1 -qi m+1) ], (A.48) 

fc =0 

i2/ = ^i? fcofcl (p (m+1) ,t)exp[^ .q( m+1) +^ 1 -qi m+1) ]. (A.49) 

fco^O 

Since R s is written as 
R s is bounded as 

\t m+1 R s \ < \t m+1 R\ < e m M. (A.51) 

Then 

\t m+1 R f \ < \t m+1 (R -R,)\< 2e m M. (A.52) 

We define 

A m+ M ( r +1 \p [m+l \t) = AUq ( r +1 \p (m+1 \t)+Rs(q ( r +1) ,p (m+1) ,t), (A.53) 

and 

B m+ M m+1 \p {m+1 \t) = R f {q {m+1 \p {m+1 \t). (A.54) 

We evaluate 

\tA m+l \ < M[ m+l \ (A.55) 

, Am+l) 

K +l {A m+ M { r +1 \p {m+l \t) - A m (q [ r +1 \p {m+1 \t)}\ < e m ^— , (A.56) 
\t m+1 B rn+1 \ < e m M^ m+1) , (A.57) 
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where 

M ( m+1 ) = M M + ( 1 + pm+i) m e m M; 

on the domain D(p m+ i). We complete the proof. 

§B Appendix; Proof of Propositions 4.2A, 4.2B 
B.l Proof of the Proposition 4.2^4 

We consider the system obtained by three times of canonical transformations. 

if« ^4 H (2) ^ iJ (3) ^ iJ (4) . 
We evaluate in how good accuracy the truncated system H$ approximates 

where 

\tM\ <M[ 4 \ \t 4 B 4 \<e 3 M^. 
By taking the difference between 

dp dB 4 
dt dq ' 

and 

we obtain 

d Ap _&b a 

dt ° dq ' 
By integrating the above inequality, we obtain 



\AP (t)\ < \AP (1)\+ Jdt^ 



< |AP (l)|+e 3 . 

If we fix 

|AP (l)|<e 3 , 

we obtain 

|APo| < e 3 
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Next we evalutate AQ X = q± — Qi- By taking the difference between 

dqi u>i f 



dt 



dA 4 dB 4 
H H -, 



and 



we obtain 



dQi Ui 1 
__ 



+ 



dA A 



d 



r <9A 4 dA 4l 



1 



+ 



+ 



d£ 4 



dt~~* L e L {oj- P y/ 2 (w-p)V2 j dPi > dpi - 

By using the mean value theorem, 

1 1 11 

(w • AP + wi ■ AP/ 



(W-P) 1 ^ ( W .P)V2 
<9A 4 <9A 4 



2 (w • P) 3 / 2 



9pi 9Pi 



9P 



9Qi 



dP^dPi 



(B.H) 
(B.12) 

(B.13) 

(B.14) 
(B.15) 



are obtained where the differentiations in the right hand side are taken at appropriate 
values between (p , z) and (P , Z). Since we are only interested in the upper bounds of 
the coefficients, we understand that the differentiations are taken at appropriate values 
between the original variables (po,z) and the truncated variables (P ,Z) without notice 
from now on. Then we get 



1 U)^ 



2 e (u-P) 3 / 2 



> • AP + u>i • APi) 



+ 



d 



d 



dP ^ dQ 1 

In the same way, as for A Pi = p\ — Pi, we obtain 



+ APi 



d 
'dP l 



3Aa dB A 



dPi dpi 



(B.16) 



d 
dt 



APi = - 



AP 



d 
dPo 



+ AQj 



•9 A „ d 
+ APi ■ 



<9Qi 



9Pi 



<9A 4 <95 4 



9Qi dqi 



(B.17) 



By using the notation z = (<Zi,Pi), Z = Pi), AZ = 2 — Z, we obtain 
A|AZ|<(^ + I)|AZ| + (^ + i)|AP | + ^, 



Under the assumption u>i ~ e, we get 



f 1 ' 



||AZ|<(I + e) |AZ| 



|APo| + F . 



By integrating the above inequality, we obtain 

|AZ| < exp [Jdt(j + e)][|AZ(l)| + J^dt\AP \ + J^dt 
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(B.18) 



(B.19) 



^B.20) 



By using (|B.10[) . and assuming 
we obtain 



lAzmi < e 3 



\AZ\ < t r exp [e{t- l)][e 3 t + e 3 ]. 
We consider the time interval as 



1 

1 < t < -s 

- - tl 3 



where 



1 



r + i 

In this interval, we obtain the evaluation as 

\AZ\ < e 2 . 

The equation for AQ = q — Q is given by 

lu> 1 



< 1. 



5** 



2 e ^.P)V^°- APo + ^- APl 

[Jn 3 . „ d A „ d ,dA 4 dB A 
+[AP ■ — + AQ 1 ■ + APa • — Ittpt + 



dQi 



OP^dPo dp 



which yields 



dt 



AQ 



<(- + h\AP \ + \\AQx\ + (1 + i)|AP x | + * 



t 



V 



t 



3 



By integrating the above inequality we obtain 

|AQ | < \ AQ (1)\ + ^ f dt\AP \ + f j\AQ!\ + J dt\AP 1 



dt- 



t 4 ' 



which yields 



AQ C 



< |AQ (l)| + i|AP | + |AZ| + e 3 

< e\ 

where we used flB.10)l . (JB.25j) and assumed that 

|AQ (1)| <e 2 . 

By obtaining (|B~To|) . (|B~25|) . (IB~29|) . we complete the proof of 



(B.21) 
(B.22) 

(B.23) 
(B.24) 

(B.25) 



^B.26) 



(B.27) 



(B.28) 



(B.29) 



(B.30) 



Next we consider in how good accuracy the perturbations of the truncated system 
Hjy approximate the perturbations of the transformed system H^ A \ We try to obtain the 
equation for A5Q$ = 5qo — SQo. We take the difference between 

1 Ll> 1 



d 

dt 6q ° 



2 e (u-pfl 2 



(u> • 5p + u>i ■ 8px) 



+(°Pa ■ -HZ - + ^i • — + opi ■ 



dp dqi 

/<- 9 . d 

+ {°Qo ■ 7T- + °Po ■ t^— + oqi 



dpi dp 



<9<?0 
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+ 5pi 



<9pi dp 



[B.31] 



and 

d . _ Iujq 1 . „ rT ^. d d d x dA A 

J t SQ " = "ST (^pF>° ' +<*•«*> + <«* ' af> + «* ■ + ■ ^ ) ^ ■ 

(B.32) 



By using the mean value theorem, we obtain 



1 CJn I 

= -0—7 • A5P o + • A5P i) 

+7 — 7 • AP o + wi • APO^o • 5P + wi • <SPi), 



(B.33) 



(5p • 77- + <tyi • 77— + 5pi • 77- ) o ( 5p o • + $Qi 



dp 



dqi 



dpi dp 



dP 



9p ttyl 

+ (AP - — + AQi 



dpi dp 



+ AP 



g 
ap 



dQi 



dP 



dP^dPo 



and 



d 
dq 



d d d ,dB 4 

opo ■ 7T— + oqi • tt— + opi ■ 



dp 



dqi 



<9Q 



(B.34) 



O r\ r\ 

= ( s Qo ■ — + Sp o ■ — + SQi ■ — + sp i 



dpi dp 

d x dB A 



dq 



dp 



dqi 



+(A5Q ■ — + A5P • — + A5Qi 



dq c 



dp c 



dpi dp 

dpi dp 



dq 



(B.35) 



where the differentiations in the right hand side are taken at the appropriate values between 
(Q, P) and (q,p), and it will not be noticed from now on. By evaluating the coefficients, 
we obtain 



< -(\A5P \ + wi|A<fPi|) + -(|AP | + u>i|APi|)(|(JPo| + wi|<SPi|) 



+(\A5P \ + \A6Q 1 \ + \A5P 1 \)- 

6 



1 



+(|AP | + |AQi| + |AP|)(|5P | + ItJQxl + |5P|)- 
+(|5Qo| + |<5Po| + |5Qi| + |<5Pi|)^ 



+(|A5Q | + |A5P | + |A^Qx| + IA5PD-. 



(B.36) 
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By using the inequalities as 



we obtain 



dt 



A5Q 



\AP \ + |AQi| + |APi| < e 2 , 
|AP | +wi|APi| < e 3 , 



< (- + 7)|A5P | + \A6P 1 \ + Ua5Q 1 \ + e -,\A6Q \ 



6 f 



t 



,3 
,3 



+e 2 |oP | + e 2 (e + 7 )|oP| + -jo^l + -\5Q 



In the same way as A5Qo, we obtain 
d 



and 



By using 
we obtain 



dt 



A5P, 



< (\A5Q \ + \A5P \ + |A(JQi| + |A5Pi|) 



+(|5Qo| + |5Po| + |5Qi| + |5Pi|)^, 



— |A5Z| < - + e |A5Z| + |AoP | + -r|AoQ 



f 1 



+e 2 (e + -)|oP | + jlSQil + e 2 (e 2 + -)|oP| + -\5Q 



|A<JQi|, |A5Pi| < |A5Z|, |5Qi|, |5Pi| < |5Z| 



^|A5Z| < (e+£)|A<fZ| + |A<SP | + ^|A<fQo| 
or t t 



+e 2 \5P \+e 2 \5Z\ + -\5Q \, 



and 



^AoQo 

dt 



< -JA5P \ + \A5Z\ + e -\A6Qo\ 



+e 2 \5P \+e 2 \5Z\ + -\5Q \, 



< (\A5Q \ + \A5P \ + \A5Z\) e - 
+(\5Q \ + \5P \ + \5Z\)^. 
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\A5Z\ < exp[ / dt{e+-)) 



{\A5Z(1)\ + J^dt\A5P \ + Jdt^\A6Q \ 
+ ^dte 2 \5P \ + J^dte 2 \5Z\ + J^dt^\5Q \} 



< t r exp[e(t - 1)] 

{lA^^I+tHA^Poll+e 3 



A5Q 



t 2 



e 2 t\\5I 



+eH\\5Z\\+t 



t 2 



}■ 



The inequality (|B.44|) yields 



\A5Q \ < exp[^tZt^] 



t, 



{|A5Qo(l)l + -||A5Po||+t||A5Z|| 



+e 2 t\\5P \\+e 2 t\\SZ\\+e 



SQ 



t 2 



} 



< {\A5Q (l)\ + -\\A5P \\+t\\A5Z\\ 

e 



+e 2 t\\5P \\+e 2 t\\5Z\\+e 



5Qt 



t 2 



and the inequality ()B.45|) yields 



|A5P | < |A5P (1)| 

A5Q 



A5Z 



If we consider the time interval 



where 



t 2 



t 2 



+ 



SPo 



t 2 



P 



r + i 



< 1, 



the growth factor of |A<5Z| is bounded as 



t r exp[e(t- l))t < -. 



34 



Then in this interval (|B.49|) . we obtain 

A5Q 



< |A5Qo(l)| + -||A5P || + l|A<yZ|| 

e 



+ e 2 \\6P \\+e 2 \\6Z\\+e 



5Qt 



t 2 



\ASP \ < |A5P (l)|+e 3 



A5Q 




A5Z 


t 2 


+ 


t 2 



+ 



t 2 



+ 



6P 



t 2 



+ 



5Z 



t 2 



+e\\5P \\+e\\6Z\\+e' 



+ e 2 


AbQo 




t 2 


SQo 




t 2 





By using (IB~53"J) to (15321) (IB~54"J) . we obtain 



A5Q C 



< |A<JQo(l)l + -|A5P (l)l + l|A5Z|| 

e 



5Q C 



-e 2 ||5Z|| +e 2 ||<5P | 



and 



||A5Z|| < -JA5Z(l)\ + -JA5Po(l)\ 



A5Q 



t 2 



5Q C 



t 2 



e||5P ||+e||5Z||. 



By using (|B.56|) to (jB.55|) . we get 
A5Q 



< \ASQ (1)\ + - e \A5P (l)\ + \\A5Z(l)\ 



6Q C 



t 2 



-e||5P ||+e||5Z|| ; 



By using (jB.57J) to ()B.56|) . we obtain 

2 



||ArtZ|| < r z \A5Q (l)\ + ^\ASP (l)\ + ~\ASZ(l)\ 



5Q C 



t 2 



e||5P ||+e||5Z||. 
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By using (IB~57)) (IB~58l) to (fB~53|) . we obtain 

||A5P || < e 3 |A5Q (l)| + |A5P (l)|+e 2 |A5Z(l)| 



5Q Q 



t 2 



+ e 3 ||,5Po||+e 3 || ( 5Z||. 



By substituting 

A5Q (1) = A5P (1) = A<5Z(1) = 
to (IB~57l) (IB~58l) JESU), we obtain 



A5Q 



5Qc 



\\A6P \\ < e 3 
\\ASZ\\ < e 2 



t 2 

5Qo 



\5P n 



t 2 

SQo 



t 2 



e 3 \\5P \\+e 3 \\5Z\ 



e||5P ||+e||5Z||. 



Next we evaluate ||<5QoA 2 || by 

d x _ 1 w 1 

= 



2 e (u>-P) 3 / 2 



(u> • SP + u>i ■ 5 Pi 



dP n 



dQi 



dPidPo 



which yields 



1 



<|5Qo(l)| + -||5P || + ||5Z||. 



By using (IB~65l) to (jESH) (lB~62j) (jB~63jl . we get 

A5Q 



(B.59) 
(B.60) 

(B.61) 
(B.62) 
(B.63) 



(B.64) 



(B.65) 



(B.66) 

(B.67) 
(B.68) 



< e 2 \5Q (l)\+e\\5P \\+e\\5Zl 

||A5P || < e 3 |5Q (l)|+e 2 ||5Po||+e 3 || ( 5Z||, 
\\A6Z\\ < e 2 |5Q (l)H-e||<5Po||+e||5Z||. 

These complete the proof of ()4.60|) . 

B.2 Proof of Proposition 4.25 

We evaluate how the errors AQ^\ AP< 4 ', A5Q {A \ A8P^ are transmitted to AQW, 
APW, A5Q^\ A5P^ by the canonical transformations 



(m+l) 



(m+l) _| 9S m 



P 



q (m) + 



dS m 

Qp(m+1) ■ 
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(m > 1) 



(B.69) 
(B.70) 



where 

\t m S m \ < e m L { ™\ (B.71) 



Lemma 



|AP (1) - AP {m) \ < ^|AQ (m) | + ^|AP (m) |, 

|AQ (1) - AQ (m) | < ||AQ (m) | + ^|AP (m) |, (B.72) 



\5P {1] -5P {m) \ < - t \5Q {m) \ + - t \5P {m) \ 



\5Q {1) -5Q {m) \ < -\5Q {m) \ + -\5P {m) \, (B.73) 

Proof of Lemma 

We can prove (1B~T2J) . (lB~73l) in the almost same way. So we prove (jB.72|) as the 

representative. We consider the difference AP = p — P. Taking the difference between 

(m) = (m+1) + J^rn (R?4) 

and 

p(m) = p(m+l) + ( g ^ 

yields 



^pH _ Ap( m +!) _|_ 



^ 9 Q( m )+^ ap(m+l)^Q(m)' ^-' D J 

where differentiations in the right hand side are taken at the appropriate values between 
{q,p) and (Q,P) according the mean value theorem. So we obtain 

| A p(m) _ A p(m+l)| < _(| A qM| + |AP (m+1) |), 

|AQ (m) - AQ (m+1) | < — (|AQ (m) | + |AP (m+1) |), (B.77) 

t m 

As for the first term in the right hand side of the above inequalities, we obtain 
| AQ (m) | < | AQ (m) - AQ (m+1) | + | AQ (m+1) | 



< — (|AQ (m) | + |AP (m+1) |) + |AQ (m+1) |, (B.78) 



which yields 



|AQ (m) | < |AQ (m+1) | + — |AP (m+1) |. (B.79) 

t 
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By using (lB~79j) to (IB~77l) . we obtain 



AQ {rn) - AQ 



(m+l) 



1 < 




AX (m+1) | 




t m 
6 '" 




1 < 




AX (m+1) | 









where 



|AX (m) | = |AQ (m) | + |AP (m) | 
In the same way as (jB.79|) . we obtain 



| A p(m)| = | A p(m+l)| + 

I II I ^ m I v 



(m+l) I 



(B.80) 
(B.81) 

(B.82) 



From (lB~79j) (lB~82i we obtain 

|AX (m) | < |AX (m+1) |. (B.83) 

We evaluate 

|APW - AP (m) | < |AP (1) - AP (2) | + |AP (2) - AP (3) | + • ■ • + (AP^ 1 ) - AP (m) | 



r m— 1 



m— 1 



~ V 1 t 2 ' 1 / 

< -|AX (m) | 

where we used (1B.80J) (IB.83J) . In the same way, we get 



|AQW _ AQ (m) | < -|AX (m) | 

These complete the proof. 
Proof of Lemma 

From the inequalities ()4.59|h we obtain 

-\AQ^\ + -\APW\ <e 2 . 



|AX (m) | 



(B.84) 



(B.85) 



(B.86) 



So 



|APW-AP( 4 )| < e 3 , 
\AQ^-AQ^\ < e 3 , 



are obtained. Therefore we obtain 



| AP (1) | < | AP W | + | AP W - AP W | < e ; 
|AZ«| < \AZ^\ + |AZ« - AZ< 4 >| < e 2 
AQ 



,(4) 



(1) 



,(4) I 



(1) 



< | AQ^ | + | AQ^_AQ^ |<e2 



t 



t 



t 



(B.87) 



(B.88) 
(B.89) 

(B.90) 
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We have proved (|4.63|) . 
Lemma 

|A5P (1) - A5P (m) | < ^(\5X (m) \\AX {m) \ + \A5X im) \), 

\A5Q {1) - A5Q {m) \ < e -(\5X {m) \\AX (m) \ + \A5X {m) \), (B.91) 

where 

\A5X (m) | = \A5Q {m) | + \A5P {m) \ . (B.92) 

Proof of Lemma 

By using ()B.79j) . we obtain 

|AQ (m) | + |AP (m+1) | < |AX (m+1) |. (B.93) 

In the same way, 

\5Q {m) \ + |5P (m+1) | < |5X (m+1) | (B.94) 
is obtained. We take the difference between 

«pW = 4- W-* ■ ^ 4- ■ (B.95) 

and 

By the mean value theorem, we get 



r) r) r) r)^ 

(SO {m) ■ +5P {m+1) - )(AO {m) ■ + AP (m+1) - V " 



<9QM Qp(m+i)^ ^ <9Q( m ) &p(™+i) y <9Q( m ) 

+(A5Q^ ■ ^ + A5P^ ■ (B.97) 

where the differentiations in the right hand side are taken at the appropriate values between 
(q,p) and (Q,P). So we obtain 

\A5P {m) - A5P {m+1) \ < (\6Q im) \ + |5P (m+1) |)(|AQ (m) | + |AP (m+1) Q — 

t m 

+(\A5Q {m) \ + |A5P( m+1 )|)^ 
< !_| ( 5x (m+1) ||AX (m+1) | 

t 

+ € Z(\A5Q^\ + \A5P^\), (B.98) 
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where (jB.93|) (|B.94|) are used. In the same way, we get 

\A5Q {rn) - A5Q {m+1) \ < ^|5X (m+1) ||AX (m+1) | 

_m 

)Hl j- I A A r>( m + 1 ) i 



+— (|A<5Q (mj | + |A5P lm+i) |). (B.99) 

From the above inequality, 

|A<5Q (m) | < \A5Q (m) - A<5Q (m+1) | + \A5Q (m) \ 

< l-|5X {m+1) ||AX {m+1) | + ^(|A5Q (m) | + |A<5P (m+1) |) 

+ |A5Q (m+1) |, (B.100) 

which yields 

|A5Q (m) | < — |5X (m+1) ||AX (m+1) | + |A5Q (m+1) | + — |A5P (m+1) |. (B.101) 
Then we get 

\A5Q (m) \ + |A5P (m+1) | < ^|<5X (m+1) ||AX (m+1) | + |A<5X (m+1) |. (B.102) 
By using f|B.1 02^ to (|B~98j) (jR99|) . we obtain 

|A<5P (m) - A5P (m+1) | < — (|<5X (m+1) ||AX (m+1) | + \A5X (m+1) \) 

m 

\A5Q {rn) - A5Q {m+1) \ < — (|5X (m+1) ||AX (m+1) | + |A5X (m+1) |) (B.103) 

t m 

By using the above inequalities, 

|A5X (m) | < |A<5X (m+1) | 

+— (\5X( m+ V\\AX( m+ ^\ + |A5X( m+1 )|), (B.104) 
t m 

is obtained. Then we get 

| A5P (1) - A5P (m) | < | A5P {1) - A<5P {2) ! + ••• + | A<5P (m - 1} - A<5P (m) | 



< -(|5X (m) ||AX (m) | + |A5X (m) |). (B.105) 



In the same way, we obtain 



|A5Q (1) - A5Q {m) \ < -(|5X (m) ||AX (m) | + |A5X (m) |). (B.106) 



These complete the proof. 
Proof End of Lemma 
From the inequalities ()4.59j) . we obtain 



AXW 



< e\ (B.107) 
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From the inequalities (jB.65|) . we obtain 



|5X (4) | < \6Q (l)\t + -\\SP \\t + \\SZ\\t. 



(B.108) 



From the inequalities ()4.60|) . we obtain 



-|A^ (4) I < e 2 |5Qi) 4) (l)l +e||5P (4) |l + 4$Z 



(4)1 



(B.109) 



By using the above three inequlities, we obtain 

\A5P^-A5P^\ < e 2 \6Q^\l)\+e\\6P^\\+e 2 \\6Z^\\ 
\A5QU-A6QW\ < e 2 |5Qi 4) (l)|+ e ||5P (4) ||+ e 2 ||5Z( 4 )|| ; 

for 1 < t < 1/e. By combining (jOUJ) (IB. HOD (IB.111I) . we obtain 



asqP 



< e 2 \8Q^(l)\+e\\SP^\\+e\\SZ^\\ 

\A6P^\\ < e 2 |5Q( 4) (l)|+e||5P (4) ||+e 2 ||5Z( 4 )| 
\A6ZW\\ < e 2 |5Q( 4) (l)|+e|| ( 5P (4) ||+e||5Z( 4 )|| 



This completes the proof of the propositionB. 



(B.110) 
(B.lll) 



(B.112) 



§C Appendix; Proof of Proposition 5.3 

The proof is essentially the same as that of the preceding proposition 4.1. So in this proof 
we estimate the generating function of the canonical transformation, and the residual parts 
which cannot be eliminated by the canonical transformation only roughly Later we omit 
m. We consider the canonical transformation induced by the generating function given by 



T l (qf\ P ^ 1 \t)=J2T k AP {l+1 \t)e^ 



(0 



(C.l) 



where T} does not depend on 

The transformed Hamiltonian is 



2 , 



m 

dt 

(J+1JU/2 



+ 



e 2 (wp('+ 1 )) 1 /2 
(0 



Ct>l • 77T 1" -fll 



(I) 



•Dt (p (m) ,t) + R 2 + E l (q[ l \ p^ ,t) + R : 



- Bl (q {l+l \p {l+l \t) + R,+ 



dt 



(C.2) 
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where R±, • • • R± are defined as 



Ri 



(u, • p^ 2 P (l+1) ) 1/2 - \~ t TTXiTTTT? • (p? - pf + \ (C3) 



e 2V ~ ^ y e 2V ~ ^ e 2 (w-p('+ 1 )) 1 /2 



R 2 = D l (jt t >,t)-D l (pl' +1 \t), (C.4) 

i? 3 = £*(<7iW) - Si(g}°,P (,+1) ,0, (C5) 

it> 4 = Bi(q {l \p il \t) - B l {q { - l+1 \p { - l+1 \t). (C.6) 

We determine the generating function T] so that the leading term depending on the slow 
angle variables can be eliminated; 

1 1 ar I ( 9 ?\p< l + 1 >,f) 



e 2 (w .p(i+i))V2 
When is decomposed as 



(0 



= -£,(qfj°,p (,+1) ,*). 



^=^e fel (p^),0e*^ 
fci^o 



(0 



7] is written as 



t,=x;^(«-p (,+i) ) 



1/2 _ 



(Wi • fci 



(0 



Since TJ is roughly estimated as 



(C.7) 



(C.8) 



(C.9) 



(CIO) 



the differences between the original variables and the transformed variables are estimated 

as 



P?-Pi' +1) = 0, 



By using these estimates, we evaluate the residual parts Ri and dTi/dt; 



where we use uj\ ~ e, 

w < 
w < 

\Ra\ < 



4 6 2 (a;-p)3/ 2lPl Pl MPl Pl j 



.2; 



t 2V 



dDi(p,t) 
dpi 



(p?- P r h , 



fi+2> 



dEi a\ n +1 \ 
"(Pi -Pi ) 



r 2l-l 



dpi 



t 2 * ' 



e 2(m-l) e « 



97j 



(C.11) 
(C12) 

(C.13) 
(C14) 

(C.15) 

(C.16) 

(C17) 
(C.18) 
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Then q$ independent (so q$ independent) part is estimated as 



i? 1 + J R 2 + J R 3 + ^-i~_ (C.19) 
while q® dependent (so qH +1 ^ dependent) part is estimated as 

2(m-l) i 

~ — T v (C20) 

The residual part 

R = R 1 + R 2 + R 3 + R4 + -^ (C21) 

is decomposed into non-oscillating part Rq, slowly oscillating part R s and fast oscillating 
part Rf] 

R = R + R s + R f , (C.22) 



where 



R = J R (p('+ 1 ),t)~-^ T , (C.23) 



R s = ^i? fcl (^ +1 ),t)exp[^ 1 . g f +1) ]~— , (C.24) 



fci^O 

2(m-l) J 

i2/ = ^^(p^^Oexp^o-^+ife!-^]-^^^. (C.25) 



feo^O 



We define 



A+i(p ( ' +1) ,t) = £>,(p (,+1) ,*) + flo, (C26) 

£m(<?f + V + V) = (C27) 

5i + i(g (,+1) ,p (,+1) ,0 = S,(flff +1 ),pf +1 \0 + «/, (C28) 

By transforming canonically, we succeeded in lowering the residual part E\ ~ e l ~ l /t l to 

§D Appendix; Proof of Propositions 5.5, 5.6 

D.l Proof of Proposition 5.5 

Lemma 



5qo 



t 



< \5q (l)\+e\5q 1 (l)\ + \\5p (l)\ + -\5p 1 (l)\, (D.l) 



e 21 " ' ■ e 



\\5 Po \\ < e 4 \5q (l)\ + e 4 \5 qi (l)\ + \5p (l)\ + e 5 \5 Pl (l)\, (D.2) 

H^ll < e 3 |^ (l)l + l^i(l)l + ^Po(l)l + l^i(l)l, (D.3) 
fell < ^\5 qo (l)\ + e 2 \5 qi (l)\+e\5p (l)\ + \5 Pl (l)\. (D.4) 
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Proof 

From the variational equations of H^ 3 ' 3 \ we obtain 



< \5q Q (l)\ + -\\5p \\ + -\\S Pl \\ + e 4 



\Spo\ < |^Po(l)|+e 4 





bO 


Sqi 


t 


1 


t 



Sqo 


+ e 4 


Spo 


+ e 4 


5qi 


+ e 4 


Spi 


t 




t 




t 




t 



Sqi 



< 1^(1)1 + -||5p || + ||5p 1 ||+e 5 



\5 Pl \ < \5 Pl (l)\+e' 



5qi 



+ 6- 



Sqo 



Spo 


+ e 2 


Sqi 


+ 6 2 


Spi 


+ e 4 


Sqo 


t 




t 




t 




t 



(D.5) 
(D.6) 
(D.7) 
(D.8) 



From these inequalities, we can deduce the inequalities of the Lemma by tedious manipu- 
lations. 

Proof End 

We assume that 



\5q (l)\ ~ |*Po(l)| ~ \Sqi(l)\ ~ \5 Pl (l)\ ~ 1, 



(D.9) 



since in the linear perturbation, the scale of the perturbation variables is arbitrary. Then 
we obtain the proposition below. 
Lemma 



Sqo 

t 


< 


1 

1?' 


(D.10) 


\$Po 


| < 


1, 


(D.11) 


Sqi 

t 


< 


1 

5 

e 


(D.12) 


\Spi 


1 < 


1. 


(D.13) 



By using the above lemma to the evolution equations of perturbation variables, we obtain 
the evaluation of proposition 5.5. 



D.2 Proof of Proposition 5.6 

Since S m ~ e 2m /t m , we obtain 



where 



|tfq (1>1) -Sq {m ' l) \ < —\8x^\ 

t 

\Sp^-Sp^\ < j\Sx^\, 



\Sx\ = \Sq\ + \Sp\. 



(D.14) 
(D.15) 

(D.16) 
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On the other hand, since T x ~ e /t and TJ does not depend on q , we obtain 

\5q M -5q {m > l) \ < e -\5y {m ' l) \, (D.17) 

|^ m,1) -^ m '°| = 0, (D.18) 

\5 P { ^ 1] -5p ( ^ l) \ < -\5y (m ' l) \, (D.19) 

t 

where 

\5y\ = \5 qi \ + \5p\. (D.20) 

By using the above inequalities, we obtain the lemma below. 
Lemma 

\5q^ - 5q^\ < € -\8y^\ + -\5qi ) m ' l) \, (D.21) 



\5p$' 1] - 5p { ™' l) \ = -\8x<> m >% (D.22) 



2 



■2 



Since 



\Sp^-Sp^ l) \ < € -\Sy^\ + j\Sqt'% (D.23) 



7^ (3 ' 3) | < -, (D.24) 
^l^ 3 ' 3) | < ^ (D.25) 
j\6x^\ < 1 (D.26) 



we obtain the proposition 5.6. 



§E Growth Index of Perturbations 

In this appendix, we calculate the growth rates of perturbations in the first model (A0^02, 
2/xi PS /i 2 ) and the second model (Xcpfcpl, ji\ ~ 112) which are presented in the beginning of 
§4, assuming that A/e is of order unity. 

The fourth order Hamiltonians of the first model and the second model are written 

by 

H = -^-P) l/2 + A + R, (E.l) 
A = gl j— ^ (M + JVcosfcQO, (E.2) 

\R\ < ^ ( R3 ) 



where 



7 = 1 k = 1 (E.4) 
M = iV = (P - 2Pi)Pi /2 , (E.5) 
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and 



7 = 2 k = 2 

M = 2N = (P — 2P 1 )P 1 , 

respectively. The coefficients of the perturbation equations 



d_ 

dt 



are given by 



/ d 2 A 



SQi 
SPi 



leu 2 



dPidQi 
d 2 A 



d 2 A \ 

2 e (u-P) 3 / 2 ' BP 2 
d 2 A 



+ 



d_ 

dt 



SQ 



dQ 2 

( 1 ^0^1 
~2 e 

V 

\u>l 1 



'w • P) 3 / 2 

• crA 

dQ^Po 
d 2 A 

+ 



dQ 1 dP 1 

d 2 A 

+ 



SQi 
SPi 



dPxdPo 



5P + R, 



SP + 



+ 



2 e (w • P) 3 / 2 <9P 2 
1 cjq^i 1 d 2 A 



+ 



5Po 



2 e (w-P) 3 / 2 9P 9Pi 
5Po(l), 



/ 

d 2 A 
dP dQ 1 

8P 1 + i?, 



SQi 



\R\ < ^(|5P | + |5Qi| + l^i|), 



(E.6) 
(E.7) 



(E.8) 



(E.9) 

(E.10) 

(E.11) 

(E.12) 
(E.13) 



r\ 1 



d 2 A 
dP 1 dQ 1 



sin /cQx H — /e 



2 (w-P) 3 / 2 



iVsinfeQi, (E.14) 



77 1_ 



d 2 A 
8P 2 



4 (w • P) 5 / 2 

1 / <9 2 M 



(M + NcoskQi) - 



( dM 



'uj.pf/2\ dPl Q Pl 



9N 

+ 7777" COS kQi 



(uj-py/ 2 yap 2 



d 2 N 
dP 2 



cos kQi 



(E.15) 



V 1 \ d 2 A 3 cj ^i Ar . _ v 1 

1 / <9M 



u; 



( dM 



2(w-P) 3 / 2 I 9Pi 9Px 



+ 7777" COS Kt^i 



9N 

+ 7777" COS Kt^i 



2(w-P) 3 / 2 yop dP 

i / d 2 M d 2 N 



(u> ■ py/ 2 \dp 1 op op^Pq 



cos kQi 



(E.16) 
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rj 1 \ 1 d 2 A 2 1 



i¥) m = - k Wpw NmskQl > (E ' 17) 

7f) 9^ = ^(^^ Wsm * gi '*(^T7^l^ sm ' ;Ql ' (E ' 18) 



and 



1( 9 2 ^ 3 ul Ar ujo ( dM ON 

77 «B2 = 77 ^M + iVcos^ ° + cosfcQj 

e f y <9P Z 4 (w ■ P) 5 / 2 (a? • P)' 3 / 2 \ <9P <9P 

+ K^( v 5Pl + 5Pj COS ^J- (K19) 

Since the correction to the growth rate of perturbations by P is of order e as long as we 
consider the finite small time range, we consider the fixed points with dropping P. 

For simplicity, we consider the case lo\ — 0. 

In the first model, around the fixed point 

Qi = | + kir 2P = P = c, (E.20) 
where k is integer, the perturbations are given by 

SQi = t ( - )fcri 5Qi(l), (E.21) 
SP 1 = t^ k+lr ^5P 1 (l) + ^5P (l))-^5P (l), (E.22) 

i 1/2 

5Qo = 5g (i)-^^p (i)(t-i) 

-1(^-1)^(1), (E.23) 

5P = <5P (1), (E.24) 

where F 1 is given by 

Fi = vtf[-L (E.25) 

The lines Qi = n/2 + /c7r and Pi = c/2 are the heteroclinic orbits. Since the definition of 
the Bardeen parameter ( contains the prefactor 1/t, when Ti is larger than 1, £ grows in 
proportion to t ri_1 . Around the elliptic fixed point 

Qi = kix 6Pi = P = c, (E.26) 

the evolution of perturbations is oscillatory. 

In the second model, the fixed points are given by 

Qi = \* 2P 1 = P = c, (E.27) 
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where for odd k hyperbolic, for even k elliptic. Around this hyperbolic fixed point the 
perturbations are given by 

fp x ) = Cl exp[r 2 (l-1)] ( V2 



c 
2 



1 



1 \ / 



+c 2 exp[-r 2 (l--)] ^2 + 5P o(l) 1 , (E.28) 



t 



2° ) V2 



where 



and 



where 



d = ^g 1 (l)-^^P 1 (l) + ^^P (l), (E.29) 
c 2 = I5g 1 (l) + ^V 1 (1)-^V (1), (E.30) 



1 1 ^ 2 

5Qo = 5Q (l)---^5P (l)(t-l) 

-¥eT^ 6 *M 1 -b + R > (K31) 

5P = SP (l) (E.32) 



|P| < ^\\SPi\\. (E.33) 



The growth rate T 2 is given by 

r 2 - - 7 -p. (E.34) 

The contribution of the rj dependent part to the Bardeen parameter ( is proportional to 
/(*) = 

/(t) = iexp[r 2 (l-i)] (E.35) 

which increases for t < T 2 and 

W = ^ exp(r2_1) ' (E ' 36) 

In case k>i 7^ 0, as the time proceeds the term originating from the unperturbed part 
becomes dominant while the perturbation parts decay as 1/t 7 , therefore the fixed points 
disappear. 
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